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Abstract By coupling observations of type la supernovae with results obtained from the best available numerical mod- 
els we constrain the Hubble constant, independently of any external calibrators. We find an absolute lower limit of H > 
50 km s^ 1 Mpc -1 . In addition, we construct a Hubble diagram with UVOIR light curves of 12 type la supernovae located in 
the Hubble flow, and when adopting the most likely values (obtained from 1-D and 3-D deflagration simulations) of the amount 
of 56 Ni produced in a typical event, we find values of H a > 66±8 and > 78±9 km s _1 Mpc -1 , respectively. Our result may 
be difficult to reconcile with recent discussions in the literature as it seems that an Einstein-de Sitter universe requires H a — 
46 km s _1 Mpc _1 in order to fit the temperature power spectrum of the cosmic microwave background and maintain the age 
constraints of the oldest stars. 
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1. Introduction 

Due to their high intrinsic luminosity and apparent uniformity, 
type la supernovae (hereafter referred to as SNe la) have be- 
come an important distance indicator in modern cosmology. 
Today, they are utilized to determine the value of the Hubble 
constant (H ) and measure the changes of the past cosmic 
expansion. For this reason a number of substantial observing 
campaigns have recentl y been co nducted for SNe la at nearby 
redshifts (see Leibundgut 2000, for a list). As a result there 
is now a considerable number of events available with superb 
temporal and photometric coverage. However, there has been 
little effort made to use these high quality data sets to link ob- 
servations with the physics of SNe la in a systematic way. The 
purpose of this article is to combine results obtained from theo- 
retical models with modern data in order to constrain the value 
afflo. 

Prior attempts to couple observations with explosion mod- 
els of SNe la in order to deter mine the value of H a include 
the p ioneer i ng inv e stigations of lArnett et alJ Jl985|). followed 
bv branch! ( Il992l). ILeibundgut & Pintol Jl992t) . iNugent et ail 
i 1995b and iHoflich & Khokhlovl dl996l) . These works gave 
promising results, constraining the Hubble constant between 
45 < H < 105 kms -1 Mpc -1 , and revealed that with 
few assumptions, SNe la used in such a manner provide an at- 
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tractive way to measure H a , while circumventing many prob- 
lems associated with the extragalact ic distance ladder (see 
b6hm-VitensJl997tlLivio et alJl997l) . 

Although several progenitor models for SNe la are dis- 
cussed, the common view favors an accreting carbon oxy- 
gen (C-O) white dwarf in a binary system, which undergoes 
thermonuclear in cineration at or near the Chandrasekhar mass 
(for reviews see iHillebrandt & NiemeverlEoOoT: lArnettlll994 
Wooslev & Weaver 1986). The energy released from burning 
to nuclear statistical equilibrium (NSE) at the density and tem- 
peratures expected in these explosions completely disrupts the 
C-O white dwarf, while the subsequent light curve is pow- 
ered by the Comptonization or "thermalization" of 7 rays pro- 
duced from the radioactive decay of 56 Ni — > 56 Co — > 56 Fe. 
The single degenerate model became favored over the double 
degenerate mod el after strong Ha emission was observed in 
SN 2002ic jHamuv et all2003t) . The Ha emission in this SN la 
is thought to come from circumstellar mater ial originating from 
the wh ite dwarf's evolved companion star faamuvetal.12 003: 
iNomotoetal 120041 

Here we use bolometric light curves of SNe la as a means to 
link observations with results obtained from models of an ex- 
ploding C-O white dwarf. Bolometric light curves constructed 
from observations provide a simple and direct route to probe 
the complicated explosion physics and radiation transport. As 
it is typically more straightforward to extract the total flux 
(hence luminosity) of a SN la from models rather than the 
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flux for individual filters, which requir e complicated multi- 
group calcula t ions jLeibundgut & Pintolll992l lEastman|ll997t 
lHoflichiri997l iLeibundeutl 12000) bolometric light curves pro- 
vide a sorely needed tool to connect observations with models. 
In addition, observed bolometric light curves are easily con- 
structed from the integration of broad-band photometry, and 
near maximum light they reflect the f raction of 7 rays ther mal- 
ized. Consequently -80% or more JSuntzefJll99a 120031) of 
the thermalized flux from the 7 rays is emitted at optical and 
near-infrared wavelengths (3000-10000 A), therefore what we 
manufacture from the observed photometry has been coined a 
UVOIR bolometric light curve. It must be noted that by neglect- 
ing the small amount of flux outside the UVOIR wavelength 
regime we introduce a systematic underestimation on the cal- 
culated values of H . However, we address this systematic er- 
ror when placing constraints on H (see below). The summed 
UVOIR flux also offers the advantage that we do not need to 
apply any K-corrections, which are necessary when just using 
individual filter passbands of SNe la located in the Hubble flow. 
Essentially, we can reduce the problem to one of energy bal- 
ances where the energy inputs from the radioactive decays and 
the losses due to 7 ray escape can be compared to the observed 
wavelength-integrated flux of the SNe la. 

In the following we utilize a set of well observed SNe la 

to demonstrate, via two methods, (see below § |4] and § 

that -under the assumption that SNe la are a product of 

the thermonuclear disruption of a Chandrasekhar-mass C- 

O white dwarf- it proves to be rather difficult to obtain a 

value of H < 50 kms" 1 Mpc -1 . Our results, along with 

recent detection of the inte grated Sachs-Wolfe (ISW) effect 
J 1 1 I . L 

(Boughn & Crittenden 2004a b) observed in the Wilkinson 
Microwave Anisotropy P robe (WMAP) data, jBennett et alJ 
l2003tlsnergel et alJl2003h cast doubts on recent discus sions in 
the literature (see Blanchardet al. 2003; Shanks 2004), which 
suggest "alternatives" to the concordance cosmological model. 
Spatially flat, matter-dominated Einstein-de Sitter models may 
produce a temperature power spectrum that can fit cosmic mi- 
crowave background (CMB) observations just as accurately as 
the best concordance model, which sports a dark energy com- 
ponent. However, Einstein-de Sitter models require very low 
values of the Hubble constant (e.g. H ~ 46 km s -1 Mpc -1 ) 
and are unable to account for the observed ISW effect. 



The structure of this paper is as follows: in §|2]we briefly 
discuss the basic details and information that we can extract 
from a UVOIR/bolometric light curve. This is followed by a 
short discussion of the data we have compiled and the method 
with which we construct our UVOIR light curves. In § 3 we dis- 
cuss different models of SNe la and the 56 Ni yields we adopt 
for a typical SN la. We then present in §E]a method to derive 
H through the combination of observations and the theoreti- 
cal 56 Ni masses calculated in explosion models. Section[5]con- 
tains a discussion of the classical way to derive H through the 
Hubble diagram of SNe la. Contrary to previous methods we 
employ here the bolometric flux. We conclude in §[6] 



2. Bolometric Lightcurves of SNe la 

2. 1. UVOIR Bolometric Lightcurves 

Here we provide a basic description of a typical UVOIR light 
curve and the physics that is thought to be driving its evolution 
and specific characteristics at different epochs. During max- 
imum light the dynamical time (i.e. time since explosion) is 
approximately equal to the diffusion time for photons trapped 
within the ejecta. This causes a reduction in the opacity, which 
then allows a larger fraction of pho tons to escape the expand- 
ing ejecta JPinto & Eastman! l200ll) . The majority of the ob- 
servable supernova flux is emitted in the optical. After maxi- 
mum light an ever increasing fraction of 7 rays escape freely 
and no longer deposit their energy in the ejecta. They are 
lost in the energy balance. In addition, due to a lack of near- 
infrared ( JH A'-passbands) data, we neglect a small contribu- 
tion of flux (no more than —5% near maximum light) in our 

I III. 
constructed UVOIR light curves (Contardo 2001). Right af- 
ter maximum light as both the 7— ray deposition rate and the 
temperature, hence opacity, decrease, it is believed that there 
is a release of 'old' photons, causing the observed luminosity 
to briefly overshoot the energy input from radioactive decay. 
After the release of stored energy the light curve declines in 
luminosity until between 20 and 40 days after maximum light, 
where an inflectio n point i s seen in most events (see Suntzeff 
Jl996h and lContardo et alJ (12000)1. After -60 days the bolo- 
metric light curve begins to follow a nearly linear decline of 
0.026±0.002 mag/day JContardo et all2000l) as the energy in- 
put from radioactive decay decreases exponentially and less 
energy is deposited by the 7 ra ys. At this time the infrared 
contribution rises to around 10% (Contardo et al. 2000). About 
a year after maximum the nebula no longer traps any 7 rays 
and the U VOIR lig ht curve is completely powered by positrons 
jMilne. The. & Leisingll999l) . 

With UVOIR light curves and accurate distances we are 
able to obtain a measure of the total luminosity and, through ap- 
plication of Arn ett's Rule, the quantity of 56 Ni produced from 
burning to NSE jArnettll982HArnett et alJl9 85). Arnett's Rule 
simply states that during the epoch of maximum light the lu- 
minosity of a SN la is equal to the instantaneous energy de- 
position rate from the radio active decays within the expand- 
ing ejecta (see alsolPinto & Eastm anll200 0albl). This rule ha s 
been utilized bylSuntzefJ Jl996h.lVacca & Leibutidgutl (fl997l) 
IContardo eUdIll2000l) . IStrolger et all J2002h . and lCandia et all 
(2003) to determine the amount of 56 Ni produced in a num- 
ber of SNe la. These efforts have revealed that the explosions 
of SNe la do indeed produce a range in the amount of 56 Ni 
synthesized from —0.1 Mq associated with the subluminous 
variety of SNe la to w 1 Mq for the most luminous ones. We 
are now in the position to use UVOIR light curves of a fair 
sample of well-observed SNe la to probe the explosion mech- 
anism. Hopefully in the near future it will be possible to place 
constraints on models as they become more sophisticated. In 
a subsequent paper we will provide a detailed analysis of the 
bolometric light curves and derived 56 Ni masses for a large 
number of well observed SNe la. 
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2.2. Observational data 

As mentioned above, there are a number of past (and present) 
dedicated monitoring programs located around the world that 
have assembled large collections of SNe la data. Programs 
which we have used here in clude: the Calan/Tololo Survey 
( Hamu v et al. I 1 19951 Il996d) the Center for Astrophysics 
jRiess et al J ll 999al llhal 12003). and the Supernovae Optical 
Infrared Survey (SOIRS) (Mamuv et all20o"fl) . 

We selected only SNe la located in the Hubble flow 
(> 3000 km s _1 ) with excellent ((7)i3Ui?7-band observations 
and that contain at least two pre-maximum observations in 
most photometric bands. Four of the SNe la compiled here 
include U -band photometry, and for those events without U- 
filter observations we added a correction (see below). At this 
stage, no corrections were made to account for contributions by 
UV-flux blueward of the atmospheric cutoff and near-infrared 
JHK-band photometry. 

Table 1 lists all the SNe la (and references) we have used 
along with information pertaining to the amount of redden- 
ing that we have adopted for each event. Values listed for 
Galactic reddening were taken from the COBE dust maps of 
ISchlegel et alJ Jl998l) . while host galaxy reddenings were pro- 
cured from a variety of literature sourc es. To be as co nsistent as 
possible we used reddenings given in Phill ips et al.l (1999) for 
all SNe la t hat coincided with ou r sample. For those events not 
included in IPhillips et all Jl 999) we adopted values from the 
literature giving preference to those calculated via the Phillips 
method. T he reddening for th e host galaxy of SN 1999dq was 
taken from lRiess et al 

Finally, in Table 1 we list the two observables that are em- 
ployed to constrain H Q . This includes the host galaxy recession 
velocity and the UVOIR bolometric flux at maximum light. 
Heliocentric velocities obtained from NED were converted to 
the CMB frame. As all SNe la listed in Table 1 are located in 
the Hubble flow, we assumed an error of 400 km s _1 for all 
velocities, in order to account for (random) peculiar motions. 
The uncertainties listed with the bolometric fluxes account for 
(1) a small measurement error, which is less than 5% and (2) 
the uncertainties associated with estimates of host galaxy ex- 
tinction. 



2.3. Construction of UVOIR Lightcurves 

We construct UVOIR light curves in the same manne r 
previously adopted b y IVacca ffi: Leibundgul ill 9961 [l997), 
IContardo etail (boOOl) and IContardol fcOOlh . The reader is re- 
ferred to these papers for a detailed discussion of this empiri- 
cal fitting method and their previously attained results; here we 
briefly summarize the main points. 

We attempt to fit SNe la photometry in a completely ob- 
jective way. Data for each filter is fitted with a ten parameter 
function. This function consists of a Gaussian, corresponding 
to the peak phase on top of a linear decline for the late time de- 
cay, an exponentially rising function for the initial rise to maxi- 
mum, and a second Gaussian for the secondary maximum in the 
VRI light curves. Fitting photometry in this manner is advan- 
tageous because a continuous representation of the light curves 



is produced without resorting to templates that may wash out 
subtleties of each filtered light curve. The ten fitted parameters 
and several other interesting quantities can b e used to e xplore 
the finer details of SNe la light curves (see IContardol 2001: 
Stritzinge rl2005l) . 

To produce a UVOIR light curve we first fit the light curve 
of each passband. Each magnitude is then converted to its 
corresponding flux at the effective wavelength and a redden- 
ing correction is applied. The flux for each filter at a given 
epoch is then integrated over wavelengths to get the total flux. 
Note, corrections are employed to account for overlaps and 
gaps between passbands . We also i nclu ded a compensation in a 
manner similar to lContardo et alJ d2000l) for those SNe la that 
have no [/-band phot ometry. IContardo et alJ used a correction 
based on SN 1994D jRichmondetal.lll995llPatatetalJll99d: 
Mei kle et alJI 199(4 ISmith et alJl2000l) 7 however, this event had 
an unusually blue color at maximum and corrections based on 
it tend to overestimate the fraction of flux associated with the 
[/-band photometry. We instead employed a correction derived 
from SN 1992A JSuntzefd lT996). which is the only well ob- 
served normal SN la with no host galaxy reddening. We esti- 
mate an additional 2% error is incurred on each UVOIR light 
curve that has our [/-band correction. 



3. 56 Ni Yields From Explosion Models 

A key ingredient for the methods presented below (see §|4]and 
§ |5j is the amount of 56 Ni produced in a typical SN la explo- 
sion. Both methods depend on th e total energy rad iated by the 
supernova to establish its distance. Contardol j200ll) showed for 
a small sample of SNe la that up to a factor of 10 difference 
in the yield of 56 Ni can exist between individual events. An 
absolute upper limit for the amount of 56 Ni synthesized is the 
Chandrasekhar mass (~1.4Mq), when the star becomes unsta- 
ble and either collapses or explodes. However, due to the pres- 
ence of intermediate mass elements (IMEs) observed in spec- 
tra taken near maximum light, we know that the white dwarf 
is not completely burned to 56 Ni. A lower limit is provided by 
the subluminous events. Although only a few of these SNe la 
have been observed in detail (due to selection effects) three 

well observed events indicate ^0.10 Mm of 56 Ni is synthe- 

i II. 
sized (Stritzinger 2005). 

To obtain a more quantitative value we turn our attention 

to recent nucleosynthesis calculations perfor med at the Max - 

Planck-Institut fur Astrophysik (MPA) ( Travag jioet al.lE"0 04). 

which are based on 3-D Eulerian hydrodynamical simulations 

dReinecke et alJ Eo02albl) of an exploding white dwarf, that 

burns via a purely turbulent deflagration flame. 1 Based on their 

highest resolution 3-D simulation (i.e. model b30_3d_768), 

which consisted of 30 "burning" bubbles and a grid size of 768 3 

for 1 octant of a sphere. iTravaglio et alJ found the total yield of 

1 Note that in 3-D deflagration simulations, once the initial condi- 
tions are set (i.e. T, p, and chemical composition) the only parameter 
that may be adjusted is the manner in which the flame is ignited. Thus 
the amount of material burned is determined by the adopted sub-grid 
model and the fluid motions on the resolved scales I Reinecke et al. 
2002a). Unlike 1-D simulations it is impossible to fine tune the amount 
of material burned at a given density. 
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56 Ni to be 0.42 Mq. However, they found that as the number of 
ignition spots is increased, more explosion energy is liberated, 
which may lead to a larger yield of 56 Ni. The number of these 
ignition spots is strictly dependent on the grid resolution of 
the simulation, which is limited by the computational facilities 
available. One may therefore expect a larger yield of 56 Ni as the 
computational power and hence grid resolution is increased. In 
addition, more recent calculations that employ ignition condi- 
tions representing a foam-like structure, consisting of overlap- 
ping and individual bubbles, indicate that it may be possible 
to liberate a larger fraction of nuclear energy as one em ploys 
different ignition conditions (Ron ke & Hillebrandi2004l) . 

We considered a range of results produced by other 
deflagration models available in the literature, in particu- 
lar the p henome nological parametrized 1-D model - W7 
jNomoto. Thielemann. & Yokoilll984T) . Recent nucleosynthe- 
sis calculatio ns show that W7 synthesizes 0.59 Mq of 
56 Ni dlwamoto et al.l ll999). It must be noted, however, that 
1-D models compared to 3-D calculations are expected to pro- 
vide a less realistic representation of the physical processes that 
occur during thermonuclear combustions because they do not 
properly model the turbulent flame physics. Additionally, mul- 
tidimensional effects are neglected which do have an important 
influence on the flame propagation. Nevertheless, W7 is a well 
establishe d mod e l that can fit the observed spectra ra ther well 
jHarknesslll99ll iMazzali et al]ll995t iMazzalilEoOll) and has 
been used extensively over the last two decades to investigate 
SNe la explosions. 

Despite the success of deflagration models, they are cur- 
rently unable to account for the more luminous SNe la, and 
predict appreciable amounts of unburned carbon, oxygen, and 
IMEs leftover in the inner ashes of the ejecta, which has 
not yet been conclusively observed. These shortcomings were 




starts as a deflagration flame until a transition occurs, causing 
the flame to propa gate supe rsonically thus the explosion be- 
comes a detonation lHoflichl i 19951) provided a series of DDMs 
that range in 56 Ni mass between ~0.34 and 0.67 Mq. His best 
fit model (M36) for the well observed SN 1994D produces 0.60 
Mq of 56 Ni. The main difference to the pure deflagration mod- 
els is that DDMs contain an additional free parameter that de- 
scribes the local sound speed ahead of the flame. This free pa- 
rameter is not physically understood but is essential to force the 
transition from deflagration to detonation. 

Throughout the following analysis we adopt results from 
the highest resolution MPA simulation and the 1-D model W7. 
Although these two models are not meant to represent the com- 
plete range in observed luminosity for the total population of 
SNe la, they produce results that are illustrative of the ma- 
jority of observed events. Both of these models are not per- 
fect and, as results shown below indicate, the MPA model may 
not be representative of the more luminous events. We take the 
56 Ni masses of these two models to be representative of a fair 
fraction of observed SNe la. 

As previously noted in § 1, just after maximum light the 
observed luminosity is expected to be larger than the radioac- 



tive luminosity, as the ejecta becomes optically thin and allows 
the release of stored UVOIR photons. At this epoch the pho- 
tosphere rapidly recedes into the ejecta, revealing deeper lay- 
ers of the pro genitor allowing more spectral lines to radiate. 
Branch ( 1992) (see his Table 1) conducted a survey of the best 
numerical models at the time and found that the models which 
adequately treat the time dependent nature of the opacity near 
maximum light predict a (the ratio of energy radiated at the 
surface to the instantaneous energy production by the radioac- 
tive decays) to be slightly larger than unity. He concluded that 
a=1.2±0.2 was the most applicable value and noted that the 
value of a appeared to be independent of the rise time. 

The parameter a=1.2±0.2 is nothing more than a cor- 
rection factor that is applied to the measured luminos- 
ity derived from the model 56 Ni masses. In the follow- 
ing, we take this parameter into account in our discus- 
sion of the values of H determined from the models. The 
56 Ni masses of 0.42 and 0.59 Mq correspond to an en- 
ergy release after 19± 3 days (the typical rise time of SNe la 



dContardo et all 120001) ) of (8.40±1.26)xl0 42 erg s" 1 and 
(1.18±0.18)x 10 43 erg s"\ respectively. If we combine the en- 
ergy production with a=1.2±0.2, the luminosity is increased 
to LOliO^xlO^ergs- 1 and 1.42±0.32x 10 43 erg s" 1 . We 
note, however, that radiatio n transpo rt calculations based on 
two MPA 3-D simulations (Blinnikov, private comm.) give the 
same value of a as calculated by Arnett's more simple analyti- 
cal models (i.e. a=l). 

Finally, we note that the value of a may vary for different 
SNe la depending on the amount and distribution of radioac- 
tive isotopes, i.e. the opacity, in the ejecta. If a=l occurs be- 
fore bolometric maximum, one would expect a smaller amount 
of stored photons. Therefore the luminosity would tend not to 
overshoot the instantaneous energy deposition rate. However, if 
a= 1 occurs after bolometric maximum, the light curve should 
in principle, be broader and flatter. 



4. H from model 5b Ni masses 

In this section we derive a first analytic expression to constrain 
H directly from model 56 Ni masses. This expression com- 
bines the UVOIR peak brightness of SNe la with explosion 
models via Arnett's Rule ( Arnett 1982). 

4. 1. Connecting H and the Model Luminosities 

First, we develop an analytic equation which uses a simple 
argument that allows one to connect H with the amount of 
56 Ni produced in a SN la explosion. This relation relies on the 
fact that H is defined as the ratio of the local expansion veloc- 
ity to the luminosity distance, which in turn is obtained from 
the inverse square law for the ratio of luminosity and the ob- 
served brightness. Therefore since the luminosity of a SN la 
depends on the amount of 56 Ni, we can use the explosion mod- 
els as our guide to the absolute luminosity. Combining this with 
both the measured brightness and recession velocity (or red- 
shift) of any particular event, we can derive a value for H . 
The first expression to constrain H therefore combines three 
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elements: (1) Hubble's law of local cosmic expansion, (2) the 
distance luminosity relation, and (3) Arnett's Rule. 
We combine Hubble's law, which is defined by 



H Q = 



cz 



with the distance luminosity relation given by 



(1) 



(2) 



where F corresponds to the UVOIR flux obtained from the 
bolometric light curves, L is the luminosity of a. fiducial SN la, 
and Dl is the luminosity distance. We obtain 



(3) 



At maximum light the luminosity produced by the radioac- 



tive ob Ni can be expressed as 



aE m (t} 



(4) 



where Em is the energy input from the decay of 56 Ni, evalu- 
ated at the time of bolometric maximum (rise time tg), and a 
accounts for any deviations from Arnett's Rule (where a — I). 
An expression for Em can be found in Nadvozhin ( 1994) 



Ni^Co 



!\i I T"~~~" - 1 ) — Qco 

A Co 



Qc D e 



-Acot 



M Ni = e(t R )M m , (5) 



where Ani and Aco are e-folding decay times of 8.8 and 1 1 1.3 
days for 56 Ni and 56 Co, respectively, and Q^i and Qc corre- 
spond to the mean energy released per decay of 1.75 and 3.73 
MeV. For 1 Mq of 56 Ni, equation l|5} turns into 



E m (l Mq) = 6.45 x io 43 e - tfi / 8 - 8 + 1.45 x 10 4 V 



•ta/U1.3 



(6) 



a normal SN la 
a peak magnitude 



iR'ess et alJ Jl999hl) found, for 
(e.g. Ami5(B) = l.l mag) with 
My = —19.45, a rise ti me to B maximum of ~19.5 
days. IContardo et alJ d2000t) found the bolometric rise time 
to be within one day of the £>-band for nearly all SNe la in 
their sample. Throughout this work we assume a bolometric 
rise time of 19=1=3 days. The adopted uncertainty should be 
adequate to account for intrinsic differences between the rise 
times of different SNe la. Using this rise time and assuming 
a=1.0 (Arnett's Rule) we can combine equations @ and 
and obtain the simple relation that gives for 1 Mq of 56 Ni a 
total luminosity at maximum light of 

L max (lM ) = (2.0 ±0.3) x lO^ergs^Mg 1 , (7) 

where the error corresponds to the 3 day uncertainty in the 
adopted bolometric rise time. 



Substituting equation (0} into equation Q we can relate the 
luminosity to the mass of 56 Ni via Em (equation|5}, and then 
if we include the factors that directly equate the luminosity 
with the 56 Ni mass, we obtain our final relation to calculate 
the Hubble constant 



H = 



( \TtF ho1 

I max 



AnF: 



hoi 



\aEm(tR) 



AnF 



hoi 



ae(tR)Mm 
(8) 



With equation (jSJl only two observables - the bolometric 
flux (-F^ax) at maximum and the redshift (z) - are required to 
determine the value of H . 

Uncertainties come from the rise time, which determines 
the the peak luminosity, the uncertainty in a, which depends 
on the radiation escape from the explosion, and of course the 
amount of 56 Ni synthesized in the explosion. Finally we note 
that our 56 Ni mass is "error free" in the sense that the adopted 
value(s) for this parameter (hence fiducial luminosities) are 
completely model dependent. 

We now have an analytic form for H which is directly con- 
nected to the 56 Ni produced in the explosion. The other param- 
eters have to do with the radiation transport: the ratio of en- 
ergy release to energy input and the time between explosion 
and maximum luminosity. 

4.2. Results 

In Fig. 1 we present results obtained using equation l|8} for all 
SNe la listed in Table 1 . For every supernova we show the de- 
rived H , assuming that its observed brightness would corre- 
spond to a given nickel mass (in steps of 0.1 Mq). The inverse 
square -root dependence of H D on the nickel mass is clearly vis- 
ible. The '1-er' error bars that accompany each point account 
for a recession velocity error of 400 km s^ 1 , an error associ- 
ated with the reddening correction, a measurement error of the 
flux (< 5%), a 3 day error for the assumed bolometric rise time, 
and a 2% error for those events that have a J7-band correction. 
Note that the most dominant error is the uncertainty associated 
with the redshift due to peculiar velocities. 

It is evident from this figure that for a given mass of 
56 Ni there exists a range of possible values of H . This is what 
we expect owing to the fact that there are known intrinsic dif- 
ferences between SNe la. Hence, if we (erroneously) assume a 
single 56 Ni mass for all observed SNe la, we obtain a range of 
H values as in Fig. 1. 

Two SNe la (SN 1992bo and SN 1993H) are clearly situ- 
ated below the rest of the obje cts and are both know n to be red 
events with Ami 5 (B) =1.69 faamuv et alJll996ct) . Both also 
show evidence of weak Ti features in their spectra JPhi llips, 
private comm.). Because the models we have adopted in this 
work were designed for normal SNe la and these two events 
are subluminous in nature, we exclude them in the following 
discussion; however, they are included in the Hubble diagram 
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(see § |5}. 2 Given the prediction of an explosion model, we 
can now read off the preferred value of the Hubble constant. 
Naively, one would like to take the mean of the distribution 
of the supernovae for a given 56 Ni mass and d erive a Hu bble 
constant. However, the natural scatter of 56 Ni (|Contardo et alJ 
l2000t IBowers et alJll997t ICappellaro et alJll997fc prevents us 
from doing this because only one value is correct for a given 
event. The best we can do now is to derive a lower limit for H 
by associating the nickel masses with the faintest supernovae 
and hence obtain an underestimate of the Hubble constant. 

By choosing a 56 Ni mass of 1 Mq and associating it with 
faint SNe la we clearly reach a lower limit for H ■ Note that 
the highest 56 Ni mas s derived from a SN la with a Cepheid 
distance is 0.84 Mq (Strol ger et alJl20 02). The solid horizon- 
tal line in Fig. 1 indicates that evidently more than 1 Mq of 
56 Ni must be produced in a normal SNe la explosion in order to 
obtain H < 50 km s _1 Mpc -1 . Note that the two dotted ver- 
tical lines indicate our adopted 56 Ni masses with a=1.0. For 
both cases an H of more than 50 km s -1 Mpc -1 is favored. 
Only the two faint SNe la fall well below this value. 

To obtain an absolute lower limit on H , we present 
in Fig. 2 the least luminous normal event in our sample - 
SN 1999aa - with 1-er and 3-cr confidence levels. This SN la 
produces values of H that are 9% below the mean determined 
from all 10 normal events. For our adopted 56 Ni masses with 
a=l we obtain a lower limit on H Q with SN 1999aa (3-cr away 
from the calculated values) to be > 40 kms -1 Mpc -1 . The 
dashed line in Fig. 2 illustrates the effect of increasing a by 
20%. We see that this gives us an additional systematic uncer- 
tainty that would decrease H by 9%. The change incurred on 
the Hubble constant lies directly on top of the quoted lower la 
confidence level. Finally, we note from equation l|8) that by ne- 
glecting ~10% of the flux emitted outside of the optical, we 
systematically underestimate H by 5%. 

We find that with white dwarfs as progenitors of SNe la it 
is very difficult to obtain a value of H Q < 50 km s -1 Mpc -1 . 
With 1 Mq of 56 Ni, one could expect H —60 km s -1 Mpc -1 . 
Observations and models for the most luminous events indicate 
that no more than 1 Mq of 56 Ni is produced. With our adopted 
56 Ni masses (0.42 Mq < 56 Ni < 0.59 Mq) we find from the 
1-cr error bars in Fig. 2 the Hubble constant to be constrain 
between 70±6 < H D < 8 3 ±7 km s -1 Mpc" 1 . 

The problem can, of course, also be inverted to derive a 
possible range of 56 Ni mass given a value of H . This will be 
interesting to constrain the 56 Ni mass for models, should H be 
known to high accuracy. For H w 70 km s -1 Mpc -1 we find 
from Fig. 1 a range in the amount of 56 Ni produced in a SNe la 
explosion to be 0.5 Mq < 56 Ni < 1.0 Mq. 

5. H through the Hubble diagram of SNe la 

With this method we determine H from the Hubble diagram 
in a manner similar to wha t has b een previously presented 
by iTammann & Leibundg utl Jl990l) and iLeibundgut & Pintol 



2 On a uniform distance scale SN 1992bo and SN 1993H are ~40% 
less luminous than the other objects and hence produce correspond- 
ingly less 56 Ni lStritzingeil2 005). 



J 19921) (see also |Sandage & Tammami 1 19931: lHamuv et alJ 
1996c; Phillip s et alJll999t IParodi et alJl2000l for similar ap- 



plications). We note that this method is similar to the previ- 
ous method, however, here H Q is calculated in a more tradi- 
tional manner. An analytic expression to constrain H D from our 
Hubble diagram is trivial to derive from the distance luminos- 
ity relation. Solving equation (|2ji for F^°l x and then taking the 
logarithm of both sides, we obtain 



log(^l) = log 



= log ( 



L,. 



\4nDl 



(9) 



Substituting cz / H for D l and rewriting the right hand side of 
equation we obtain 



log(F™J = -21og(cz)+log(L r . 



l0g(47T)+2l0g(J? ) 

(10) 



There is a linear relation between \og(F r °l x ) and log(cz) as 
can be seen in Fig. 3. A linear regression to the data yields a 
slope of 2.01 ±0.25, which is fully consistent with the expected 
slope of 2 for a linear local expansion derived in equation (II 0i . 
From Fig. 3 it is also obvious that the two faint objects are 
outliers. If they are ignored, the fit sharpens up to a slope of 
1.97 ± 0.10. With a fixed slope to this linear expansion value 
of 2 we derive the y-intercept, which corresponds to 



log(L, 



■log(47r)-21og(iJ ). (11) 



Solving equation il Q for H , we arrive at our final expression 
for the Hubble constant 



H = 



-in 



(12) 



The Hubble constant is now simply calculated by plugging 
in the y-intercept, b, derived from the linear regression of 
the Hubble diagram and a fiducial luminosity defined by our 
adopted models. 

5. 1. Results 

In Fig. 3 we present our Hubble diagram containing all SNe la 
listed in Table 1 . Error bars for all events account for both a 
redshift uncertainty of 400 km s -1 and the uncertainties asso- 
ciated with host galaxy reddening. A weighted least-squares fit 
to the Hubble diagram (for all 12 SNe la), with a fixed slope 
of 2, yields b =3.292±0.047. Plugging this into equation O 
along with our adopted 56 Ni masses of 0.42 Mq and 0.59 Mq 
we find Ho to be > 85±7 and > 72±6 km s -1 Mpc -1 (l- 
<t error) respectively. Accounting for a=1.2±0.2 we obtain 
lower values of the Hubble constant to be H > 78±9 and > 
66±8 km s -1 Mpc -1 , respectively. 

As an upper limit in the amount of 56 Ni synthesized in the 
most luminous SNe la explosion is expected to be ks 1 Mq, 
we can use the corresponding luminosity as a guide to obtain 
a lower limit on the Hubble constant through the Hubble dia- 
gram. Thus, for 1 Mq of 56 Ni and a—1 we obtain a value of 
H > 55±5 kms" Mpc -1 . As in the previous method, we 
are underestimating H by —5%, due to flux outside the optical 
window. 
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6. Discussion 

Under the main assumptions: (1) that the progenitors of SNe la 
are C-0 white dwarfs that explode at or near the Chandrasekhar 
mass, and (2) that the amount of 56 Ni synthesized to first or- 
der determines the peak luminosity, we are able to use results 
obtained from state-of-the-art numerical simulations of explo- 
sion models to uniquely define the bolometric peak luminos- 
ity, and in concert with photometric observations constrain H . 
The attractiveness of this approach stems from the ability to 
bypass assumptions that are typically made when one attempts 
to determine H , i.e. the extragalactic distance ladder and its 
accumulation of error from rung to rung. We stress that our fit- 
ting method does not add any corrections to the data. In other 
words we d o not nor malize the flux to any decli ne rate rel ation 
(e.g. A mi. JPhillips et alJll999t) . MLCS jRiess et al.lll996l) or 
stretch JPerlmutter et all 19971) V 

Uncertainties that marginally affect (in decreasing order of 
importance) our results include: the abundances of peak Fe 
group metals - stable and radioactive - produced in the ex- 
plosion models, the redshift peculiar velocities of each SN la, 
the total absorption, the assumed rise time to bolometric maxi- 
mum, the exact nature of a, which may slightly vary from SN 
to SN depending on the exact nature of the opacity and ioniza- 
tion structure of the expanding ejecta, and the amount of flux 
that we neglect outside of the optical window. 

It is still unclear what parameters effect the amount of 
56 Ni produced in a SN la explosion. Obvious candidates are the 
initial conditions prior to explosion. These include the metal- 
licity of the C-O white dwarf, the central density and the ig- 
nition mechanism. If there is a considerable fraction of al- 
pha elements such as 22 Ne within the progenitor one would 
expect more stable isotopes such as 58 Ni and 54 < 56 Fe to be 
prod uced from burning to N SE, thus reducing the 56 Ni yield 
(e.g. iBrachwitz etldl EoOO). A higher central density on the 
other hand would lead to a more robust explosion and hence 
an increased amount of 56 Ni. As discussed earlier the explo- 
sion mechanism is still uncertain and different deflagration 
and detonation scenarios produce different amounts of peak 
Fe group elements. Nevertheless with a larger 56 Ni mass we 
obtain smaller values of H . Errors attributed to the adopted re- 
cession velocity and reddenings produce scatter in our Hubble 
diagram but have very little effect on our calculations of H via 
equations (jSJi or d 1 21 . The ±3 day departure from our adopted 
rise has no more than a 10% effect on our calculations. A 
change in a by ±20% can affect H by 9%. Finally, we reit- 
erate from equation (jSJl that by neglecting ~10% of the flux 
emitted outside of the optical, we are underestimating H Q by at 
least 5% or correspondingly more if more flux is unaccounted 
in our method. 

In § 15.11 we determined a rather high value 
(85 kms _1 Mpc _1 ) for the Hubble constant when using 
results from the MPA model. This indicates that the amount of 
56 Ni produced in these 3-D deflagration simulations currently 
are on the low side. And indeed a large sample of SNe la show 
that the average distribution of 56 Ni mass is slightly higher 
(~ 0.6 Mq) llStritzingeri2005l) . This suggests that their models 



may need more fine tuning in order to produce a larger amount 
of 56 Ni, and hence match observations mo re accurately. 

There have been many attempts since pre- 
sented his Hubble diagram to exp l oit SN e la to determine 
H . We refer the reader to Branch (1998) for a detailed re- 
view of previous works that attempt to calculate H D based 
on SNe la. He concluded, from methods based on physical 
considerations similar to the methods presented in this work 
and methods which utilize SNe la that have been indepen- 
dently calibrated by Cepheids, a range in the Hubble con- 
stant of 54 < H < 67 km s _1 Mpc -1 , with a "consen- 
sus" on H = 60±10 kms" 1 M pc " 1 . More r ecent investi- 
gations of ISuntzeff et alJ (^999) and Ijhal J2002t) give values 
of H = 64 kms" 1 Mpc" 1 . Finally. JFreedman et alJ l200lt 
ISpergel et all 120031: lAltavilla et ail 12004) have all measured 
slightly larger values of H « 70 km s" 1 Mpc _1 with 10% ac- 
curacy. 

Another method independent of the extragalactic distance 
ladder which combines X-ray imaging of galactic clusters with 
the Sunyaev-Zerdovic h effect (SZE) has been recently used 
to place limits on H„ JMvers et alJll997t iMason et alJEoOll: 
Ijones et alll2003t iReese et alJl2002l) . These works have pro- 
vided the detailed study of 41 clusters giving distances which 
yield an a veraged valu e of H « 61±3 km s _1 Mpc _1 (for a 
review see Rees el2003l) . 

We find from both methods presented here that the Hubble 
constant must be > 50 km s _1 Mpc _1 in order to be compati- 
ble with current supernova models. In addition, we stress that 
this lower limit is based on the assumption that 1 Mq is an up- 
per limit on the amount of 56 Ni produced in a SN la explosion, 
and not from our adopted models. This result, along with other 
methods to measure H D using SNe la calibrated with Cepheids, 
SZE/X-ray distances and evidence of the ISW effect, strongly 
suggest that we do not live in a matter dominated universe with- 
out some form of cosmological constant or similar agent. 
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Table 1 . Well-observed SNe la in the Hubble Flow 



SN 


Filters 


Ref. 


E(B-V)^ 


E(B-V) host 


VCMB b 

(kms^ 1 ) 


p6oI 
A max 

(erg s _1 cm~ 2 ) 




SN1992bc 


BVRI 


1 


0.022 


0.000 


5870 


(1.565±0.124)xl0~ 


n 


SN1992bo 


BVRI 


1 


0.027 


0.000 


5151 


(9.106±0.967)xl0 _ 


12 


SN1993H 


BVRI 


1 


0.060 


0.050 


7112 


(4.640±0.613)xl0" 


12 


SN1995E 


BVRI 


2 


0.027 


0.740 


3478 


(5.726±0.622)xl0 _ 


1 1 


SN1995ac 


BVRI 


2 


0.042 


0.080 


14651 


(3.425±0.477)xl0~ 


12 


OIN 1 77JUU 


-D V HI 


9 
z 


n AOs 


n i sn 

U. lju 




to s/t?+n s9iwin _ 

\A. JH-Z.ZIZU. JL 1 } A ±U 


11 


SN1996bo 


BVRI 


2 


0.078 


0.280 


4857 


(2.382±0.237)xl0~ 


11 


SN1999aa 


UBVRI 


3, 4,5 


0.040 


0.000 


4546 


(2.333±0.446)xl0~ 


11 


SN1999aw 


BVRI 


5,6 


0.032 


0.000 


11754 


(3.525±0.700)xl0" 


12 


SN1999dq 


UBVRI 


3 


0.024 


0.139 


4029 


(3.871±0.730)xl0~ 


11 


SN1999ee 


UBVRI 


7 


0.020 


0.280 


3169 


(5.781±0.881)xl0~ 


11 


SN1999gp 


UBVRI 


3, 8 


0.056 


0.070 


7783 


(9.270±1.082)xl0~ 


12 



a Taken from lSchlegel et aljfl998T) dust maps. 

b Heliocentric velocities from NED transformed to the cosmic microwave background frame. To ac- 
count for peculiar velocities we assume throughout this work an error of 400 km s" 1 for all CMB dis- 
tances. 
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Figure 1 Hubble constant as a function of 56 Ni mass for 12 
SNe la in the Hubble flow for a= 1 . Error bars represent all 
uncertainties discussed in section RTTl The dotted lines indicate 
the two 56 Ni masses from our adopted models. The solid verti- 
cal line indicates an upper limit on the production of 56 Ni ex- 
pected from the thermonuclear explosion of a C-O white dwarf, 
and the solid horizontal line corresponds to the lower limit we 
derive for H . 



Figure 3 Hubble diagram of 12 SNe la in the Hubble flow. The 
negative logarithm of the UVOIR flux at maximum light is plot- 
ted against the logarithm of the recession velocity in the CMB 
frame. Error bars account for a peculiar velocity of 400 km s ~~ 1 
and uncertainties associated with host galaxy reddening. 
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Figure 2 To highlight the absolute lower limit we derive on H 
we plot results obtained with SN 1999aa (the least luminous of 
our 10 normal SNe la) with l-a (dark shading) and 3-cr (light 
shading) confidence levels. Vertical solid lines indicate both the 
adopted 56 Ni masses and the 1 Mq upper limit expected from 
an exploding C-O white dwarf. The dashed line illustrates the 
effect if a is increased by 20% (i.e. a=1.2). 



